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Abstract. In a recent paper (Usoskin et al., 2002a), we have reconstructed the concentration of the cosmogenic 
^'^Be isotope in ice cores from the measured sunspot numbers by using physical models for ^"Be production in 
the Earth's atmosphere, cosmic ray transport in the heliosphere, and evolution of the Sun's open magnetic flux. 
Here we take the opposite route: starting from the ^''Be concentration measured in ice cores from Antarctica and 
Greenland, we invert the models in order to reconstruct the 11-year averaged sunspot numbers since 850 AD. 
The inversion method is validated by comparing the reconstructed sunspot numbers with the directly observed 
sunspot record since 1610. The reconstructed sunspot record exhibits a prominent period of about 600 years, in 
agreement with earlier observations based on cosmogenic isotopes. Also, there is evidence for the century scale 
Gleissberg cycle and a number of shorter quasi-periodicities whose periods seem to fluctuate in the millennium 
time scale. This invalidates the earlier extrapolation of multi-harmonic representation of sunspot activity over 
extended time intervals. 



1. Introduction 

Sunspot numbers (SN) form the most common index of 
solar activity reflecting the varying strength of the hy- 
dromagnetic dynamo process which generates the solar 
magnetic field. 

There are two approaches to reconstruct SN for times 
prior to regular direct observations. The first approach is 
based on mathematical extrapolation using the statisti- 
cal properties of the SN record observed during the last 
300 years. Such models provide, e.g., a single (11.1-year) 
carrier frequency or a multi-harmonic representation of 
the measured SN, which is then extrapolated backward in 
time HSchove 19551 |Nagovitsyn 1997[[Rjgozo et al. 2001| ). 
Some models also use fragmentary data from naked-eye 
sunspot observations or sightings of aurorae. The disad- 
vantage of this approach is that it is not a reconstruction 
based upon a description of physical processes but rather 
a prediction based on extrapolation. Clearly such models 
cannot include periods exceeding the time span of obser- 
vations upon which the extrapolation is based. Hence, the 
pre- or post-diction becomes increasingly unreliable with 
extrapolation time and its accuracy is hard to estimate. 
The second approach is based on measured archival 
of SN such as cosmogenic isotopes 
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Fig. 1. Raw and smoothed ^''Be data. The lower curves 
(left axis) give the raw (thin curve) and the 1-2-1 filtered 
8-year- averaged data from Antarctica IjBard et al. 1997)l . 
The upper curves (right axis) show the raw and 
the 11-year smoothed yearly data from Greenland 
HBeer et al. 1990|l . 



proxies ot blN such as cosmogenic isotopes <j.g. 
ILingenfelter 1963| IBeer et al. 19901 lO'Brien et al. 1991 
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IDamon fc Sonett 1991llBeer 2000|l . While this approach is 
based upon real measurements, the quantitative relation- 
ship between SN and the measured cosmogenic isotope 
concentration is usually described in an ad-hoc fashion by 
a simple inverse relation. 

In this paper we follow an improved version of the 
second approach and reconstruct the SN for the last mil- 
lenium using a physical model with a realistic relation 
between the concentration of the cosmogenic ^"Be isotope 
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Fig. 2. The differential response function of ^°Be to the 
two main species of galactic cosmic rays, protons and 
a— particles, for different values of the heliospheric mod- 
ulation strength, <&. 

in polar ice and the SN. The ^"^Be concentration corre- 
sponds well, with a response time of 1-2 years, to the 
flux of primary cosmic rays in the vicinity of the Earth 
HBard at al. 19971 IBeer 2000|) . We use two sets of ^^Be 
data shown in Fig. ^ One is the annual series of ^"Be 
concentration in Greenland ice (Dye-3 site, 65.15 N 43.82 
W) for the years 1424-1985 HBeer et al. 1990|l . The other 
series gives the ^"Be concentration with roughly 8-year 
averaging in Antarctic ice (South Pole) for the years 850- 
1900 l|Raisbeck ct al. 1990llBard et al. 1997|l . 

2. Model 



In a recent paper ljUsoskin et al. 2002a|l . we have devel- 
oped a model which connects the sunspot number N with 
the cosmic ray flux (CR) through the source term, S, the 
open solar magnetic flux, Fg, and the modulation strength, 
through a sequence of steps: 



Be 



(1) 



Steps (1) and (2) are based upon a model for the 
open solar magnetic flux IjSolanki et al. 2007)1 . steps (3) 
and (4) on a ID model of heliospheric transport of 
CR l|Usoskin et al. 2002b|l . For step (5) we have earher 
H Usoskin et al. 2002a|l assumed that the ^°Be production 
and the related concentration in polar ice is proportional 
to the flux of CR with an energy of about 2 GeV. Here we 
calculate step (5) using a more realistic production rate 
of ^°Be by CR in the Earth's atmosphere. The local ^"Be 
production rate, R, is given by 



/■OO 

R = J X(P) ■ Y{P)dP, 



(2) 



where Pc, X{P), Y{P), and P are the local rigidity cutoff, 
the differential CR spectrum near the Earth, the differen- 
tial yield function of ^'^Be production, and the primary 
CR particle's rigidity, respectively. 
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Fig. 3. Production rate of ^°Be in the atmosphere as a 
function of the heliospheric modulation strength, and 
the local rigidity cutoff, Pc- 



Heavier species (mostly a— particles) compose a part of 
all CR, about 6% in particle number or about 25% in nu- 
cleon number according to recent high-precision measure- 
ments dg^^rr^ IBoenzio et al. 19991 lAlcaraz et al. 2000al 
IAlcara^r^^2000b|l . Moreover, since heavier nuclei are 
less modulated in the heliosphere than protons, their 
contribution is larger at lower energies (about 1-2 
GeV/nucleon) which are most effective for ^°Be pro- 
duction in the atmosphere. Therefore, we have taken 
a— particles into account in the present study. Species 
heavier than a— particles can be neglected because of their 
low abundance in CR). 

Using differential CR (protons and a— particles) 
spectra obtained from our heliospheric model 
HUsoskin et al. 20 02^)1 . and the differential yield 
function of ^°Be production in the atmosphere 
dWebber fc Higbie 20031 ), we have calculated the i°Be dif- 
ferential response function (Fig. which is the integrand 
appearing in Eq. (jSJ. One can see that the contribution 
of heavier nuclei to ^"^Be production becomes significant 
at energies lower than 2 GeV/nucleon, while the overall 
production rate R is mostly determined by protons. The 
mean energy for ^'^Be production move in the course of a 
solar activity cycle, varying between 1.5 GeV (<i>=300 MV 
corresponding to recent solar cycle minima) and 4 GeV 
($=1000 MV corresponding to solar cycle maxima), which 
is somewhat larger than the value of 2 GeV used earlier 
HMcCracken fc McDonald 20011 lUsoskin et al. 2002ajl . 
This results from the fact that the magnitude of the solar 
cycle modulation of CR greatly increases towards lower 
energies, which decreases the effective energy responsible 
for the ^°Be variations with respect to the mean energy. 
The total ^°Be production rate is shown in Fig. as 
a function of the heliospheric modulation strength, <I>, 
and Pc- The resulting ^"Be concentration in polar ice is 



I. G. Usoskin, K. Mursula, S. Solanki, M. Schiissler and K. Alanko: Solar activity: Millennium scale 



3 




10000 



Fig. 4. Scatter plot of 11-year averaged (iV^) vs. (N)'^ 
based upon the yearly group sunspot numbers for 1700- 
2000. The regression line has a slope of 1.32 ± 0.01. 




Fig. 5. Actual and restored yearly group sunspot numbers 
for 1700-1985 as a test of the consistency of the inversion 
model. 



assumed to be proportional to the ^''Be production rate, 
i.e., we ignore the details of the atmospheric transport 
and deposition of ^°Be (.Masarik fc Beer 1999|l since we 
are interested in much longer time scales. 

We now invert the model by going through the se- 
quence of Eq. ^ in reverse order: 



N 



(1') 



S 



(2') 



Fo 



(3') 



(4') 



CR 



(5') 



^Be 



(3) 



Steps (5') and (4') are taken in one stride so that from 
a measured ^"Be concentration in polar ice at a given 
location (given Pc), we determine the corresponding he- 
liospheric modulation strength, $, using the relationship 
shown in Fig. 01 Steps (l')-(3') are inversions of the corre- 
sponding steps (l)-(3) in the direct model. In step (3'), a 
power-law relation between $ and the open magnetic flux, 
Fo, is used IjUsoskin et al. 2002a|l : 



Fo = 0.023 



(4) 



where $ and Fo are given in MV and 10^^ Wb, respec- 
tively. 

In step (2') we invert the following equation which con- 
nects the evolution of the open magnetic flux Fo with the 
source function S IjSolanki et al. 2000)l . 



dFo _ Fo 
dt r' 



(5) 



where r = 4yr represents the characteristic decay time. 
This equation includes differentiation, and when employ- 
ing it to compute S from Fo, its inversion is stable only for 
smooth (noiseless) data series, for which the time deriva- 
tive does not fluctuate strongly from one data point to the 
next. Real ^°Be data are noisy and do not fulfill these con- 
ditions (see Fig.^. Therefore, a reconstruction of yearly 
SN values from the yearly sampled Dye-3 Greenland ^"Be 
record is extremely noisy. However, averaging the yearly 
Greenland data by an 11-year running mean prior to re- 
construction greatly reduces the noise and increases the 
stability of the procedure, although at the cost of los- 
ing information on variations shorter than a solar cycle. 
Similarly, a direct reconstruction using the 8-year-sampled 
South Pole record is also somewhat noisy. We therefore 



smooth the Antarctic data set using a 1-2-1 filter, which 
roughly corresponds to an 11-year smoothing. In what fol- 
lows we only consider the (effectively 11-year) averaged 
data series, so that Eq. (0) takes the form 



at T 



(6) 



where the angular brackets indicate temporal averages. 
From the time sequence of averaged measured ^°Be data 
we determine through steps (3')-(5') the open flux values, 
{Fo). The corresponding values of the average source func- 
tion, (S), are obtained by replacing the time derivative in 
Eq. 10 by a simple first-order finite difference. 

The source function, S, is the following function of the 
sunspot number iV 



S{N) = a (24.35 + 22iV - 0.0617V^) 



(7) 



where a = 1.95 • 10" Wb/yr IjSolanki et al. 2000|l . In 
order to obtain the average sunspot number, (N), we 
have to to invert the time-averaged Eq.Q and to ex- 
press {N'^) in terms of (iV). Unfortunately, there is no 
unique solution to this problem. A regression analysis 
based upon the 11-year- averaged group sunspot number 
(GSN) record Poyt fc Schatten 1998| ) since 1700 leads to 
a factor 1.32 ±0.01 (Fig.^. Therefore, we have used the 
relation (N^) = 1.32(iV)2 in the inversion of Eq.(CI). 

In order to check the consistency of the full inver- 
sion procedure we have performed the following test. 
Using the direct model (Eq. QJ, we calculated, simi- 
lar to IjUsoskin et al. 2002a|l but with the more realis- 
tic ^°Be production described above, the expected ^°Be 
concentration from the actual yearly group sunspot data 
( |Hoyt fc Schatten 1998| ). The calculated ^°Be concentra- 
tion was then used as an input for the inversion (Eq. O , 
from which the corresponding restored group sunspot 
numbers were determined. The result is shown in Fig. 
Since the restored SN series is very close to the actual 
record (the RMS errors are 10.4 for the yearly and 2.4 for 
the 11-year averaged series), we conclude that the inver- 
sion procedure is consistent. We emphasize that a nearly 
perfect reconstruction such as that in Fig. is only possi- 
ble for almost ideal data, with little noise, good sampling. 
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Fig. 6. 11-year averaged sunspot numbers for the last 400 
years: group sunspot numbers (GSN), Wolf sunspot num- 
bers (WSN), SN reconstructed from Greenland ^"Be and 
SN reconstructed from 1-2-1 filtered Antarctic ^°Be data. 
The shaded area roughly indicates the range of averaged 
SN for which the model is less reliable. 

etc. The measured -'^^Be data do not fulfill these conditions 
nearly as well as the artificial record ^°Be constructed from 
GSN, as discussed above. 

3. Reconstruction of sunspot numbers 

As a first step, we have reconstructed the sunspot activity 
for the last 400 years from ^°Be data using the model 
described above. Since direct sunspot observations exist 
for this period, we can compare the reconstructed and 
measured SN. The SN reconstructed from the two ^"Be 
series are shown in Fig. together with the measured 
Wolf and group sunspot numbers. The agreement between 
the reconstructed SN series based on the 11-year averaged 
Greenland ^°Be and the measured SN is quite good for 
the period 1700-1985, with an RMS discrepancy of 11 
(the cross-correlation is r = 0.78^!^:?^) for GSN and 12.5 
(r = 0.74:to°l) for the Wolf sunspot series. The RMS 
deviation is 10 between the 1-2-1 filtered 8-year averaged 
GSN and SN reconstructed from the Antarctic ^°Be series 
for 1700-1900 {r = 0.761^5:52) ■ 

We note that the overall agreement between the re- 
constructed and measured SN is somewhat better for the 
group sunspot series than for the Wolf series. For exam- 
ple, in 1770-1800 and 1830-1880, the reconstructed SN 
is very close to the GSN while the WSN is systemati- 
cally higher. The only exception is the period 1700-1750 
when the reconstructed SN shows good agreement with 
the Wolf series while GSN is systematically lower (Fig. ^ . 
This might imply that the group sunspot number values 
are somewhat too small during the first half of the 18*'' 
century, in general agreement with the results of Letfus 
(2000). We note that the two reconstructed SN series lie 
closer together than the GSN and WSN series in the com- 
mon time interval, i.e., between 1700 and 1900. 

In summary, we conclude that our method recon- 
structs the measured sunspot numbers reasonably well 
for the last 400 years. The only exception is the 
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Fig. 7. Sunspot activity (11-year averaged) during the 
pre-instrumental era: earlier extrapolations by Schove 
(1955 - S55), Nagovitsyn (1997 - N97) and Rigozo et al. 
(2001 - ROl), and our SN reconstructions from 1-2-1 fil- 
tered 8-year Antarctic (An) and from 11-year averaged 
Greenland (Gr) ^°Be series. The shaded area indicates 
the range of averaged SN for which the applicability of 
the model is limited. Horizontal bars with arrows denote 
the known grand minima and maxima: Maunder minimum 
(Mm), Sporer minimum (Sm), Wolf minimum (Wm), Oort 
minimum (Oni) and the Medieval maximum (MM). 

Maunder minimum. Whereas the actual sunspot number 
is close to zero ||Eddy 1976| IRibes fc Nesme-Ribes 19931 
Usoskin et al. 2000|l during this period, the reconstruc- 
tion returns values between 10 and 25. This is due to the 
fact that the model is not reliable during extended peri- 
ods of very low solar activity, as was already pointed out 
l|Usoskin et al. 2002?ji|l . On the other hand, the SN is still 
well reconstructed during the Dalton minimum circa 1800. 
Accordingly, we estimate that the reconstructed sunspot 
values are unreliable (overestimated) when SN< 20. 

In the next step, we have reconstructed the SN in the 
millennium time scale (Fig. EJ) using both the Greenland 
(since 1424) and the Antarctic (since 850) ^°Be time se- 
ries. Great minima as deduced from proxies like ^^C, 
^"Be, etc., also appear in our reconstruction, as does the 
Medieval maximum IjStuiver fc B raziun as 1989)l . We note 
that, because of the above-mentioned limitations of our 
model, the reconstructed SN, although significantly re- 
duced during great minima, never reach zero. While the 
SN reconstructions based upon the Antarctic and the 
Greenland data are rather consistent with each other af- 
ter 1600 (Fig. El, they exhibit notable differences in 1480- 
1600. The Antarctic series shows an extended minimum 
while the Greenland series exhibits pronounced maxima 
in 1500 and in 1560. Since this interval corresponds to the 
Little Ice Age, climatic differences between Greenland and 
Antarctica may well be the source. 

Fig. [7| also shows that for the pre-instrumental 
era before 1600 our reconstruction differs significantly 
from earlier results obtained by extrapolating multihar- 
monic representations of the measured SN IjSchove 19551 
[Nagovitsyn 1997] [Rigozo et al. 2001] - henceforth referred 
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Fig. 8. Autocorrelation function of the measured SN 
(group, GSN, and Wolf, WSN) and the reconstructed 
(from Greenland and Antarctic ^°Be data) sunspot series. 
Left panel: autocorrelation function for the period 1700- 
2000 (1700-1900 for the Antarctic series; 1700-1985 for 
the Greenland series). Right panel: autocorrelation func- 
tion for the period after 850 for the Antarctic series and 
after 1424 for the Greenland series. Error bars represent 
the 95% confidence interval. 
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Fig. 9. Power spectra of the sunspot number as recon- 
structed from the Antarctic (850-1900) and Greenland 
(1425-1980) ^"^Be series and of the group sunspot num- 
ber (GSN) series (1700-2000). The right panel is a zoom 
of the hatched rectangle in the left panel. Thin lines in- 
dicate the corresponding 95 % confidence intervals. Bars 
in the upper panels represent the basic frequencies used 
by S55, N97 and ROl in their multi-harmonic models of 
sunspot activity. 



to as S55, N97 and ROl, respectively). For example, the 
S55 and N97 series exhibit high maxima of activity around 
1100, 1400 and 1550-1600, which are either not present 
or significantly lower in our reconstruction. On the other 
hand, ROl predicts a high maximum around 1050 when 
the Oort minimum, visible in all other series, is expected. 
Throughout the whole period covered by Fig. [7| our re- 
construction shows much smaller variations of the SN level 
than the earlier series. 

4. Periodic features 

As discussed above, our reconstruction recovers the great 
minima and maxima of solar activity, which have been 
qualitatively deduced earlier. We now analyze the periodic 
features of the reconstructed SN series. Since we deal with 
either 11-year averaged or 1-2-1 filtered 8-year sampled 
series, we will discuss only periods longer than 40 years. 

First we compare the periodicities in the reconstructed 
series with those in the measured SN after 1700. In Fig.|Sl 
(left panel) we plot the autocorrelation function (ACF) for 
this period for the two reconstructed and two measured 
SN series. The pattern is the same in all the four series, 
the main feature being a period of about 100 years, the 
Gleissberg secular cycle. One can see that all the curves 
are rather close to each other, within the 95% confidence 
interval, giving further support to the validity of our re- 
construction method. The right panel of Fig. |H1 shows the 
ACF for the last millennium for the reconstructions from 
the Greenland and Antarctic ^°Be data. Two dominant 
periodicities are present in the reconstructed data. A cy- 
cle with an approximate period of 130 years is indicated 
by the consecutive peaks at 130, 260, 400, 530, . . .years 
in both SN series. 

Moreover, a cycle with a period of about 600 years 
visible as a deep minimum with the lag of about 300 years 



in both series and as a maximum at 600-650 years in the 
Antarctic series. 

Fig. ini presents the power spectra of the two 
reconstructed SN series and the GSN series. The 
95 % confidence levels of the power spectrum 
l| Jenkins fc: Watts 1969|> are shown by thin lines. Note 
that the period after the Maunder minimum is too 
short to have a reliable power spectrum for the 8-year 
sampled Antarctic series, so that we have plotted the 
power spectrum using the whole data set. The power 
spectra for the GSN and for the Antarctic data are in 
a fair agreement in the period range from 40 to 100 
years, revealing similar main periods: 52, 65-75 and 100 
years in GSN series and 45-50, 65 and 90-95 years in 
the Antarctic series. The peak at 130 years, however, is 
not found in the GSN record. The power spectrum of 
the Greenland series (1425-1980) shows a very strong 
peak at 60 years which is absent in the other two series. 
Other significant peaks (50, 80 and 100-115 years) 
correspond roughly to similar peaks in both GSN and 
Antarctic power spectra. Periods longer than half of 
the analyzed interval cannot be reliably determined as 
shown by thin lines in Fig. The power spectrum of the 
Antarctic series also exhibits a peak at about 220 years 
probably related to the so-called de Vries or Suess cycle, 
l|Suess 19801 [Wagner et al. 2001)) and a p rominent peak 
close to 600 years ( |Sonett fc Finney 1990| ). The latter is 
in agreement with the autocorrelation function shown in 
Fig. IHl At the top of Fig. 1^1 the periods used by previous 
investigators to extrapolate the SN are marked. Not all 
of them agree with peaks in the power spectra of GSN 
and our reconstructions. We also note that periods of 
different quasi-periodicities (such as, e.g., the Gleissberg 
cycle) seem to be fluctuating on the long-time scale. This 
invalidates the earlier extrapolation of multi-harmonic 



I. G. Usoskin, K. Mursula, S. Solanki, M. Schiissler and K. Alanko: Solar activity: Millennium scale 
90 



180 




1000 1200 1400 1600 1800 2000 

Years 

Fig. 10. Variation of the parameters of the geomagnetic 
dipole during 850-2000 AD. Migration of the dipole axis 
(top panel, in geographical coordinates) and change of the 
virtual dipole moment M (lower panel). 

representation of sunspot activity over extended time 
intervals. 



5. Geomagnetic variations 

The above results were obtained assuming a local ^"Be 
deposition and a constant rigidity cutoff of Pc = 0.5 
GV in Eq. which roughly corresponds to the atmo- 
spheric cutoff. The geomagnetic field is known to change 
in time both in strength and orientation, which affects 
the ^°Be production rate on time scales of centuries and 
longer | |Baumgartner et al. 1998| ). Fig. 1101 shows the vari- 
ation of the geomagnetic field for the period 850-2000 
AD (see also Fig. 7 of [Hongre et al. 1998| ). The param- 
eters of the geomagnetic field were taken from Hongre 
et al. (1998) for 850-1700 AD, from Bloxham & Jackson 
(1992) for 1700-1900, and from the IGRF (International 
Geomagnetic Reference Field) model for 1900-2000 AD. 
Using these data, we calculated the expected variation of 
the vertical geomagnetic cutoff, Pc, at the two -^"Be sites 
using St0rmer's equation l|Elsasser et al. 1956jl : 



Pc(xM ■ cos'^A, 



(8) 



where M is the virtual dipole moment and A is the local 
geomagnetic latitude. The calculated value of Pc remains 
negligibly small (below 0.025 GV) for the South Pole dur- 
ing the whole interval studied here. The geomagnetic lat- 
itude of the South Pole site remains above 80°, i.e. in- 
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Fig. 11. Global production rate of ^°Be in the atmosphere 
as a function of the virtual geomagnetic dipole moment, 
M, and the modulation strength, $. 



side the polar cap where geomagnetic field lines are open. 
Therefore, only the atmospheric cutoff was applied for the 
South Pole. The value of Pc for the Dye-3 Greenland site 
remains below 0.5 GV during the period of available data 
(since 1424) and, therefore, can be neglected as well. At 
earlier times, however, the influence of the geomagnetic 
field becomes significant, e.g., the value of Pc exceeds 1 
GV for the Dye-3 site before the 14*'' century. 

We have also calculated the global ^"Be production 
in the entire atmosphere as a function of the magnetic 
dipole moment and the modulation strength (Fig. see 
also Fig. 2 in IBeer 2000|l . As an extreme case, we have 
performed a reconstruction of the SN assuming global 
mixing of ^"Be before deposition, similar to -'^''C. The re- 
sult is shown in Fig. 1121 and is compared with the re- 
construction based upon local ^"Be deposition and the 
GSN record. Since the geomagnetic dipole was stronger 
in earlier times leading to a more effective shielding of 
the Earth, the same amount of globally distributed ^°Be 
would require a stronger CR flux and, correspondingly, 
reduced sunspot activity. Therefore, the SN level corre- 
sponding to the global ^"Be model is systematically be- 
low the local model and also the GSN level. This differ- 
ence becomes increasingly larger for earlier times and the 
sunspot level goes systematically below 10 before 1600 in 
the global ^°Be model. This implies that the local deposi- 
tion model is closer to reality than global mixing, although 
the truth probably lies in between. Furthermore, the local 
^"Be model gives an upper limit to the level of sunspot 
activity since in the presence of atmospheric mixing the 
geomagnetic effects would require reduced activity levels 
at earlier times in order to retain consistency between the 
^"Be and GSN records. 
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Fig. 12. 11-year- averaged sunspot group simspot number 
series (GSN) and the reconstructed SN for local and global 
^"Be deposition are represenetd by grey, solid, and dotted 
curves, respectively. 

6. Conclusions 

The sunspot number record is not only the most widely 
used proxy of solar activity, it is also one of the longest 
running scientific data series. Nonetheless, for many pur- 
poses it is still too short. For example, both the statistical 
investigation of the solar dynamo and studies of a possible 
influences of solar activity on climate would greatly profit 
from a longer time series. 

Since direct observations are sparse and not very reli- 
able before the era of telescopic observations we have used 
an indirect proxy, the concentration of the cosmogenic 
^°Be isotope in the ice sheets in Greenland and Antarctica. 
Unlike previous investigators, however, we have employed 
physical relations to convert the ^"^Be concentration into 
averaged sunspot numbers. A comparison of directly ob- 
served SN with the reconstructed values for the last 400 
years confirms the reliability of the latter. To enhance the 
reliability over longer times we have also considered the 
possible influence of geomagnetic fleld variations. Taking 
the available data with a time resolution of the length of 
the solar cycle we have reconstructed the SN back to 850 
AD. 

We have also performed a detailed analysis of period- 
icities in the long-term solar activity record. Both applied 
methods, autocorrelation function and power spectrum, 
reveal that a number of strong periods in the directly mea- 
sured sunspot series are also present in the SN record re- 
constructed for the same period of time. Some differences 
between records are found, but they are generally smaller 
than changes in dominant frequencies as different lengths 
of the time series are considered. In the higher frequency 
range, periods of about 50 and 65 years are visible in both 
measured and reconstructed series. The secular Gleissberg 
cycle, which is roughly 100 years in the GSN series, is split 
into the dominant 120-130-year and subdominant 90-95- 
year cycles in the reconstructed series, the latter being 
more prominent during the last 3-4 centuries. There is 
also evidence for an independent periodicity of about 220 
years. This may be related to the de Vries (also called 
Suess) cycle, earlier known in, e.g., -^'^C isotope data. A 



strong cycle with a period of roughly 600 year is apparent 
in the reconstructed series based upon the Antarctic data. 
Both the GSN record and the the SN series reconstructed 
from the Greenland ^''Be data are too short to exhibit this 
period, which is possibly related to the 650-year periodic- 
ity in the ^''C data IjDamon fc Sonett 199l)l . 

In conclusion, we have presented here a new recon- 
struction of solar activity on the millennium time scale 
based upon a description of the related physical processes. 
The results will be the subject of further analysis. 
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